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COLOR SLIDE SUPPLEMENT: IN- 
STALLMENT ONE 

For those of you who have been wondering when to ex- 
pect the supplement you have paid good money for, know 
that I expect to send out the first 12 slides with Amy #7, and 
13 slides subsequently. 

All of these have been created by John Dewey Jones, who 
is bound to be at least slightly disappointed if his slide sub- 
scription yields nothing but copies of his own stuff. Any- 
one have good color slides they’d like to submit? 


THE FARSIDE — THE MANDELBROT 
SET ON THE HP 28C 

Only a true follower of Richard Nelson (1) would ever 
think of evaluating the Mandelbrot set on a pocket calcula- 
tor. Richard was the somewhat non-linear guru who led a 
happy band into regions of the Hewlett Packard calculators 
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where H-P had not particularly intended them to tread. Send- 
ing these machines to unlikely levels of activity was encour- 
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aged. I plead some justification for the following. Antici- 
pating fast floating point support and color on the Macin- 
tosh, I started to look at Turbo Pascal for the Mac. Any 
self-respecting Mac program is called up by clicking the 
mouse on an icon. An icon is a 32 x 32 pixel digital image, 
and the screen of the HP 28c happens to be 32 pixels high. 

While one could obviously just draw an icon, my distinct 
cast of purism demanded an authentic calculation. 

The other excuse is that the 28c is a particularly versatile 
machine and can deal with a complex number as easily as a 
real one. Although the Mandelbrot algorithm is scarcely 
complicated to program, on the 28c the central loop becomes 
short and simple. Procedure MAND takes a complex num- 
ber associated with a pixel of the image onto the machine’s 
unlimited parameter stack, written in the form (Re,Im). It 
passes back on the stack the modulus or size of the iteration 


of z4 27 +0, and increments the counter. 
MAND 
<< SQ C + 'count' 1 STO+ DUP ABS >> 
This procedure is repeatedly called by LOOP, which 
tests the size and the count. 
The language has affinities to FORTH and to LISP. 
One could do a more detailed version on the printer, but 
since even this miniature takes a couple of hours it does not 
seem worthwhile. The HP 41c is probably an order of mag- 
nitude slower. 
Henry Casson (5047) 


(1) This is not really a footnote; members of PPC append 
their serial number to their name like apprentice robots. 
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PLOTTING THE MANDELBROT SET: AN 
INCREASE IN EFFICIENCY 

It has been been conjectured, and it seems in practice to be 
true, that any area totally surrounded by pixels of one count 
value is internally of that value. 

I make use of this to speed up generation of plots. Track 
round the edges of a rectangular area, computing the iteration 
count for each pixel. If all these are the same, the area can 
be immediately filled. If not, the procedure calls itself recur- 
sively with the longest axis halved, and this repeats until the 
rectangles become of minimum size. 

What follows is sketchy Pascaloid pseudocode. 


Procedure ProcessPix 
takes in the X and Y coordinates of the pixel, computes 
the real numbers for its position on the complex plane, 


246 to 


and returns the iteration count for the usual z + z 
the array PixelArray[1..Width,1..Height]. 

Prefill PixelArray with a flag value defined as ‘Unpro- 
cessed’. This must be outside the range of allowable iter- 
ation values returned by ProcessPixel. 

Procedure BinaryTree(Left, Top,Right,Bottom); 
begin 
if Right-Left > 1 { criterion to terminate recursion } 
then begin 

if PixelArray[Left+1,Top-1] = Unprocessed 

{ the rectangle has not been filled in a previous 
call } 

then begin 

AllSame := TRUE; { Boolean flag } 
if PixelArray[Left, Top] = Unprocessed 
then ProcessPix(Left, Top); 
For | := Left+1 to Right do { top of rectangle } 
begin 
if PixelArray[l, Top] = Unprocessed 
then ProcessPix(I, Top); 
if PixelArray[I, Top] + TestValue 
then AllSame := FALSE; 
end 
For | := Top-1 downto Bottom do { Right side } 
begin 
if PixelArray[Right,!] = Unprocessed 
then ProcessPix(Right,|); 
if PixelArray[Right,l] + TestValue 
then AllSame := FALSE; 
end 
For | := Right-1 downto Left do 
begin 
if PixelArray[|,Bottom] = Unprocessed 
then ProcessPix(!,Bottom); 
if PixelArray[I,Bottom] + TestValue 
then AllSame := FALSE; 
end 
For | :=Bottom+1 to Top do 
begin 
if PixelArray[Left,l] = Unprocessed 
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then ProcessPix(Left,|); 
if PixelArray[Left,] + TestValue 
then AllSame := FALSE; 
end 
if AllSame { We have been all round the edges} 
then FillRectangle(Left, Top,Right, Bottom, 
TestValue) 
else if Right-Left = 2 { A single line of pixels } 
then for | := Top-1 downto Bottom+1 do 
ProcessPixel(Left+1,]1); 
if Right-Left < Top-Bottom 
then begin 
BinaryTree(Left, Top, Right, 
(Top+Bottom) Div 2); 
Binary Tree(Left, 
(Top+Bottom) Div 2,Right,Bottom); 
end 
else begin 
BinaryTree(Left, Top,(Right+Left) Div 2, 
Bottom); 
BinaryTree((Right+Left) Div 2, Top,Right, 
Bottom); 
end 
end 
end 
end 
Henry Casson 
3843 §.W. 58 Drive; Portland, OR 97221 
The picture on the following page was also sent in by Dr. 
Casson. He notes that it was plotted as 769 x 1025 pixels 
using a 2 pen Hewlett-Packard 7470A driven by a Gridcase 


II, (4.77 MHz 8086) + 8087. The calculation took 8 hours 


using Turbo Pascal with the main loop coded inline, and the 
plot took 1.5 hours. The program is available on 3.5" IBM 
disks. 


MANDELBROT MACHINE DESIGN 
From Ulrich Schmidt (July 26, 1987): 
Dear Rollo, 

Here is a block diagram of my Mandelbrot machine (see 
page 4). The floating-point chips (adder, multiplier and reg- 
ister caches) are clocked at 16 MHz, while the control logic 
(sequencer, address generator) ticks along at 8 MHz. One mi- 
crocode word is therefore sufficient to supply the control 
signals for two floating-point clock cycles. 

The writable control store is made up of static CMOS 
RAMs organized 4Kx4. These chips are from Cypress; they 
have a 45 ns access time and separate pins for input and out- 
put — that helped to save a lot of glue logic. Another eight 
of them were used for the dual-port buffer. The arbitration 
logic is very simple: the first half of a 125 ns clock cycle is 
dedicated to the host whereas the other half is reserved for use 
by the Mandelbrot machine. 

(Continued on page 4) 
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The host machine I use is a 68010 box from Eltec (a Ger- 
man VMEbus manufacturer) running under OS9/68K. The 
system has 2 MB RAM on the main board and another 2 MB 
used as a frame buffer for a graphics sub-system. Hitachi’s 
fine 63484 does all the hard work, and a 1024x768 color 
monitor with look-up table is used to view the results. The 
VMEbus is a very nice asynchronous 32 bit bus which can 
be easily interfaced to a synchronous system such as the 
Mandelbrot machine. However, any other system could be 
used as host instead (imagine: I have a 68020/68881 board 
running as a peripheral on the venerable Apple II!). 

I wire-wrapped the board — there is simply nothing that 
beats wire-wrapping in areas such as cross-talk avoidance, 
wire length minimization, and ease of use. A Motorola en- 
gineer once told me that they had a wire-wrapped 25 MHz 
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68020 board running trouble-free in their lab! Okay, let me 
admit that I made a virtue of necessity: with four of these 
power-hungry pin grid arrays (144 pins) there is not much 
choice — either use wire-wrapping in combination with a 
well-soldered power distribution system, or design a multi- 
layer board with lots of ground and Vcc planes. The latter al- 
ternative costs a fortune during development and makes sense 
only if you can sell the board in volume. I honestly doubt 
whether the “Mandelbrot market” is that large. So, this 
board is not intended to be a commercial product. Itis possi- 
ble to reproduce it from the schematics if and only if you 
have access to a very good logic analyzer (400 MHz) and a 
good scope (4 channels, 300 MHz). 

The nearest thing to a microcode assembler I have seen is 
the “Meta-Assembler” sold by STEP Engineering. As the 
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STEP assembler is a rather expensive product (around $5000) 
I decided to write my own assembler to run under OS9/68K. 
I have someone coding it for me, and it should be ready by 
the time you read this. In the meantime, I hand-assemble 
my microcode algorithms, which is almost as much fun as 
toggling in computer programs from the front panel. Still, 
the result is worth it, as anybody will attest who has seen it. 

You surely want to know how your DivCon algorithm 
fares on my machine. Well, I must disappoint you: except 
for very special cases (such as all points having the same 
dwell) DivCon is significantly slower than TV scan. For ex- 
ample: TV scan computes the base picture (-2.25-1.125i, 
0.754+1.125i, 1024x768, dwell = 256) in 29 seconds, while 
DivCon needs 82 seconds. The ratio is even worse for more 
interesting areas such as the seahorse valley. The reason for 
this is that my board is essentially a vector processor, and 
vector processors are best at simple things. Anything that 
interrupts the pipelined operation degrades performance sig- 
nificantly — much more so that with a normal microproces- 
sor. You see, there are a lot of decisions to be made in Div- 
Con (Do the perimeter pixels have the same dwell? Should 
we partition horizontally or vertically? Have we reached the 
base of the recursion?). Even worse, we cannot overlap com- 
putations because we have to wait for the results in order to 
determine how to carry on. All this excess baggage is not 
there in TV scan; host machine and vector processor can 
work in parallel — while the vector processor computes the 
nth horizontal line, the host reads the n-1st line from the 
dual-port buffer and sends it to the graphics sub-system, pix- 
el by pixel. 

Where do we go from here? My board’s 16 MFLOPS 
could be upgraded to 40 MFLOPS if the new 3210/3221 
chips from Analog Devices were used. However, for now 
those chips are just paper tigers, and even if they can be pur- 
chased some years from now they will probably require half a 
power station to run. The Weitek chips are real, and they do 
work. I haven’t had the chance to try out AMD’s new 29300 
series, which looks quite interesting, too. Is anybody at a 
VLSI design centre reading this? Here is a nice little student 
project: design a chip or two which implement complex addi- 
tion and multiplication in hardware. Give the user a choice 
between single and double precision. Forget about IEEE 
compatibility (except for the format) and concentrate on 
speed instead. Use low-power CMOS (please!). The result 
would be a “Mandelbrot chip(-set)”. 

Best regards, 

Ulrich Schmidt 

An der Junkersmiihle 33/35 
5100 Aachen, West Germany 
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56. (Cover picture) Science News 131,25 (June 20, 1987). 
“The sound of a fetus breathing can be converted into a sym- 
metrized dot pattern, as shown. This particular technique, 
which depicts sound waves so that frequency information is 
emphasized, is one of several methods recently developed by 
IBM researcher Clifford Pickover for representing complex 
data. Such representations can be used to enhance subtle dif- 
ferences between waveforms and to identify patterns and 
trends. (Illustration: Pickover/IBM)” 


57. “Picture This”, Science News 131,25 (June 20, 1987) 
392-395. [Discusses Clifford Pickover’s work in depicting 
speech sounds and other complex data in novel ways.] 


“COMPLEX ANALYSIS” 
Some of you may be interested in getting a good introduc- 
tion to complex analysis. It would be hard to find a better 
one than: 

Complex Analysis (3rd Ed) 

Lars V. Ahlfors 

McGraw-Hill Book Company (1979) 

ISBN 0-07-000657-1 

The price is $40.95 — a bit steep for a book of only 331 
pages — but I found it worth the price. If any of you want 
it, I can probably get it for $37.00, postpaid: let me know. 


LETTERS 
From Mark Bolme 
Dear Rollo: 

We are putting the finishing touches on FractalMagic for 
the Mac. After writing our assembly routines, I reviewed 
yours and I had some comments that you may want to incor- 
porate into your routines. When rounding register D7, you 
typically doa 
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ADD.L #$8000,D? 
CLR.W D? 
SWAP D? 


After the ADD.L, you need to check the carry bit and per- 
form suitable operations should D? overflow. One choice is 
to 

GERALD? 

ADD.L #$8000,D? 

CLR.W D? 

SWAP D? 

BCC.S Continue 

ADDX.L D7,other data register 


When forming 2ab 
ADD.L D7,D7 


you also have to check the overflow. It is more serious here, 
because 50% of the time the result will be wrong if you 
don’t. Try ASL.L #1,D7 instead, but you still need to check 
overflow. Generally, you need to check the overflow on all 
your adding operations. 

When you form the B29 scaling from the B26, I believe 
that you simply do an ASL.L (I didn’t check too closely). If 
I am right, this will lose two bits of information in the 
number. You need to keep track of the full 64 bits of the 
product, and capture the two top bits of the bottom 32 bits 
of the product when converting into B29. Keeping all 64 
bits will also make dealing with overflows easier. I hope 
these comments are of some interest to you! 

Mark W. Bolme, President 

SINTAR SOFTWARE 

P.S. I didn’t check your code very closely. If my sugges- 
tions don’t make sense, my apologies! 


R.S. replies: 

Receiving Mark’s letter plunged me into a careful scrutiny 
of the assembler code in my EMD program. Mark is indeed 
correct about the loss of two bits of information due to the 
shift. I have recoded the dwell routine to deal with this. It is 
not, however, necessary to check the carry bit and take care 
of overflow. A careful analysis of the program reveals that 
such a condition can never occur — for rather obscure rea- 
sons. 

Because of this bug, I am no longer distributing EMD2. 
It will be corrected in the next release — EMD 3. 


From André van Meulebrouck 
Dear Amygdala: 

I wanted to write to tell you that I appreciate your publi- 
cation, I’m proud to be a subscriber, and I intend to be a sub- 
scriber for quite some time to come. Perhaps in the future I 
will submit an article to Amygdala on fractals, mathematics, 
and/or plotting fractals in LISP (MacScheme specifically, 
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which is what I write all my plotting programs in). 

In the meantime however, I have much to learn about 
fractals — especially the Mandelbrot set. I will soon be or- 
dering the book “The Beauty of Fractals”, and perhaps that 
will answer many of my questions that the Mandelbrot book 
did not answer (I think I need something more tutorial). 

Nonetheless, being extremely busy, I don’t have as much 
time to devote to fractals as I would like, and therefore I es- 
pecially would like to request more Sesame Street style arti- 
cles. A good example of what I’m suggesting is shown in 
the article from issue #1a: “Complex Numbers: the Basics”. 
I really liked that article and will be studying it more deeply 
as time permits. I would bet that articles even simpler, 
aimed at rank beginners would be greatly appreciated by your 
subscribers. 

I feel that the primary goal of Amygdala should not be to 
impress, but to humbly educate, and I suspect that many 
subscribers are subscribers precisely to get such information 
to help them get started, and to provide continual tutorial 
help thereafter. 

For example, I still do not know exactly now to plot 
Mandelbrot fractals. An article that walks a person through 
it and assumes nothing about previous knowledge of com- 
plex numbers or even plotting would be nice. 

Also, I would enjoy seeing information on random fields 
fractals, and the other sets from time. It seems a bit narrow 
to focus in on one set of the 4 major types (random fields, 
Koch curves, dragons, Mandelbrot set). 

Finally the M-fiction is a nice contribution, but not at the 
expense of tutorial articles! 


R.S. replies: 

I appreciate your comments, André. Please let me know 
if Henry Casson’s program shows you how to plot Mandel- 
brot fractals, and if there are specific areas in which you 
would benefit from more detailed treatment. We have a wide 
variety of mathematical skills represented among our readers, 
and that tends to create problems. 


COMMERCIAL PRODUCTS 
If you place an order as a result of seeing the follow- 
ing notice, please mention Amygdala with your order. 


ART MATRIX; PO Box 880; Ithaca, NY 14851 
USA. (607) 277-0959. Prints, FORTRAN program 
listings, 36 postcards $7.00, sets of 2 packs $10.00, 140 
slides $20.00. Or send for FREE information pack with 
sample postcard. Custom programming and photogra- 
phy by request. Make a bid. 


CIRCULATION 

As of August 23, 1987 Amygdala has 159 paid-up sub- 
scribers, of which 41 have the supplemental color slide sub- 
scription. 
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